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We derive a novel model-independent result for the pion susceptibility in QCD via the isovector- 
pseudoscalar vacuum polarisation. In the neighbourhood of the chiral-limit, the pion susceptibility 
can be expressed as a sum of two independent terms. The first expresses the pion-pole contribution. 
The second is identical to the vacuum chiral susceptibility, which describes the response of QCD's 
ground-state to a fluctuation in the current-quark mass. In this result one finds a straightforward 
explanation of a mismatch between extant estimates of the pion susceptibility. 
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Colour-singiet current-current correlators or, equiva- 
lently, the associated vacuum polarisations, play an im- 
portant role in QCD because they are directly related 
to observablcs. The vector vacuum polarisation, e.g., 
couples to real and virtual photons. It is thus ba- 
sic to the analysis and understanding of the process 
e+e^ -^hadrons [1, 2]. In addition, analysis of the large 
Euclidean-time behaviour of a carefully chosen correlator 
can yield a hadron's mass [3, 4]; and correlators are also 
amenable to analysis via the operator product expansion 
and are therefore fundamental in the application of QCD 
sum rules [5]. 

In the latter connection, the vacuum pseudoscalar sus- 
ceptibility (also called the pion susceptibility) plays a role 
in the sum-rules estimate of numerous meson-hadron cou- 
plings; e.g., the strong and parity- violating pion-nucleon 
couplings, g-nNN and /ttNN, respectively [6, 7, 8]. Fur- 
thermore, as will become plain herein, the pion suscepti- 
bility is as intimate a probe of QCD's vacuum structure 
as the scalar susceptibility [9] but its veracious analysis 
is more subtle, with conflicts and misconceptions being 
common [7, 8, 10, 11, 12]. 

We approach the vacuum pseudoscalar susceptibil- 
ity via the isovector-pseudoscalar vacuum polarization, 
which can be written^ 



'{P;0=N,tYZ4 



l75r'5(<j+)*ri(q;F)5(g_), (1) 



where the trace is over flavour and spinor indices; C is the 
renormalisation scale; Z4(C,A) is the Lagrangian mass- 
term renormalisation constant, which depends implicitly 



on the gauge parameter;^ and J d'^q/{2iT)'^ rep- 

resents a symmetry-preserving regularisation of the inte- 
gral, with A the regularisation mass-scale which is taken 
to infinity as the last step in a complete calculation. 

Herein we will subsequently assume isospin symmetry; 
viz., equal u- and d-quark current-masses, in considering 
the isovector-channcl. An extension to three flavours and 
the flavour-singlet channel can be pursued following the 
methods of Ref. [13]. 

In Eq. (1), S is the dressed-quark propagator and Fs 
is the fully-dressed pseudoscalar vertex, both of which 
depend on the renormalisation point. The propagator is 
obtained from QCD's gap equation; namely, 

S{p)-^ = Z2(z7-p + mi^-) + I](p), 

g'^D^cuip - q)—l^,S{q)—T^{q,p), 

(2) 

where D^y(/s) is the dressed-gluon propagator, T^{q,p) 
is the dressed-quark-gluon vertex, and m'^'" is the A- 
dependent u- and d-quark current-quark bare mass. The 
quark-gluon-vertex and quark wave-function renormali- 
sation constants, Zi,2{C,t also depend on the gauge 
parameter. 

The gap equation's solution has the form 



s{p)-^ - i-^-pA{p'-e)+B{p'-e) 



(3) 



and the mass function AKjP') = B{p'^,(^^)/A{p'^,(^) is 
renormalisation point independent. The propagator is 
obtained from Eq. (2) augmented by a renormalisation 
condition. 



^ In our Euclidean metric: {7^,7,^} 
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=^ P^ <0. 
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Physical quantities obtained from Eq. (1) arc manifestly 
invariant. 
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Since QCD is asymptotically free, the chiral limit is 
defined by 



Z2(C,A)m'^-(A) = 0, VA»C, 



(4) 



which is equivalent to requiring that the renormalisation 
point invariant current-quark mass is zero; i.e., to = 0. 
A mass-independent renormalisation scheme can then be 
implemented by fixing all renormalisation constants in 
the chiral limit [14]; namely, one solves the chiral limit 
gap equation subject to the requirement 



■ p . 



(5) 



This is implicit in the subsequent analysis. We note that 
Z2(C,A)to'^-(A) = Z4(C,A)to(C), (6) 

where rn{Q is the familiar running current-quark mass. 

The pseudoscalar vertex is determined from an inho- 
mogeneous Bethe-Salpeter equation; viz.. 



^4 [-75 T-^]tu - 



(7) 



where k is the relative- and P the total-momentum of the 
quark-antiquark pair; r, s, t, u represent colour, flavour 
and spinor indices; 



xiik;P)^S{k+)Ti{k;P)S{k^), 



(8) 



k± = k zL P/2, without loss of generality owing to 
the symmetry-preserving nature of the regularisation 
scheme; and K{q,k;P) is the fully-amputated two- 
particle-irreducible quark-antiquark scattering kernel. 

Much of the preceding material recapitulates results 
familiar from QCD's Dyson-Schwingcr equations (DSEs) 
[15, 16], which also provide the foundation for our sub- 
sequent analysis. Consider, then, that in the presence 
of a spacetime-independent pseudoscalar source, 5*5 ^ 0, 
associated with a term 



d xq{x) -15T-Ss,q{x) 



(9) 



in the action, one can define a vacuum pseudoscalar con- 
densate, whose gauge-invariant, properly renormalised 
form in QCD is 

/•A . 

a^5(s5, to; C, A) = Z^N^iv j '-j^r'Siq; sg, ™; C) • (10) 
This is analogous to the vacuum quark condensate [17] 



cr(TO;C, A) = ZiNcti- 



ir°5(q;TO;C) , 



(11) 



~ diag[l, 1], whose source-term in the QCD action is 
that associated with the current-quark mass. 

It should be emphasised that when to = 0, it is only 
foreknowledge of nonzero current-quark masses via the 



Higgs mechanism which leads one to express dynamical 
chiral symmetry breaking (DCSB) as 



(qqfc 



lim (j{m;(,A) ^ 0, 



(12) 



instead of choosing a different vacuum vector; e.g., 
(cr;s5) oc (0; 1,-1,0). Moreover, Eq. (12) defines what 
we mean by an isoscalar-scalar configuration: isovcctor- 
pseudoscalar correlations are by convention measured 
with respect to this configuration. 

These observations highlight the importance of the 
pseudoscalar susceptibility 



_d_ 



■cr5(s5,m;C, A) 



(13) 



S5=0 



Following the method of Ref . [9] , it is straightforward to 
show that 



(0 = -2 uA^ {P = 0; S5 = 0, to; C) . 



(14) 



NB. Whilst hitherto we have not specified a regularisa- 
tion procedure for the susceptibility, it can rigorously be 
defined via a Pauli-Villars procedure [9] . 

We are interested in the value of X^-' (C) in the neigh- 
bourhood of the chiral limit. Therein one may write [18] 



iri{k;0) = hj,r^El^{k;0) + 



-ri(fc;0), (15) 



where r^^{k; P) is the pion bound-state's canonically- 
normalised Bethe-Salpeter amplitude; ?'7r(C), determined 
by 



iVctrZ4 



(16) 

l^,T''xi{q;P), (17) 



is the residue of this bound-state in the inhomogeneous 
pseudoscalar vertex; and E§{k] P) is a part of the in- 
homogeneous pseudoscalar vertex which is regular as 
P^ + ml-^ 0. 

These statements will not be surprising once one re- 
calls that the solution of a linear, inhomogeneous integral 
equation is a sum; viz., the regular solution of the inho- 
mogeneous equation plus a solution of the homogeneous 
equation, here, naturally, the canonically normalised so- 
lution. In terms of this solution, the pion's Icptonic decay 
constant is expressed through 



NAtZ-? 



'qW) (18) 
h5l^.r\i{q;P). (19) 



One can turn to the axial- vector Ward-Takahashi iden- 
tity in order to determine E^{k; P); viz.. 



P^n,ik;P) 



2to(C) tTiik;P) 



S-\k+)h-f,T^ + U-f,T^S-\k_)., 



(20) 
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where rg^(fc;P) is the inhomogeneous axial- vector ver- 
tex. At P = with m ^ there is no pole contribution 
on the left-hand-side and hence Eq. (20) states 



i{OEi{k-P = Q)^B{e-m-e), 



(21) 



namely, this regular piece of the pseudoscalar vertex is 
completely determined by the scalar part of the m ^ 
quark self-energy. NB. It is straightforward to ver- 
ify Eq. (21), order-by-order, via the gap and Bethe- 
Salpeter equations using the systematic, nonperturba- 
tive, symmetry-preserving DSE truncation scheme intro- 
duced in Refs. [19, 20]. 

We insert Eq. (15) into Eq. (14) to obtain 



TABLE I: Vacuum pseudoscalar susceptibility and related 
quantities, computed using the two kernels of the Bethe- 
Salpeter equation described in connection with Eqs. (30), 
(31 and (32). Dimensioned quantities are listed in GeV, 
K ~ ~{{<iQ)1Y^^ Ei-nd /tt is the pion's chiral-limit leptonic de- 
cay constant. The entries were compiled from Refs. [9, 21]. 
NB. For quantitative comparison with some other studies 
[7, 8, 10, 12], our results for X should be multiplied by (27r)^. 



Vertex 










m 






RL [Eq. (31)] 


1 


1 

2 


0.25 


0.091 


0.0050 


1.77 


0.39 


BC [Eq. (32)] 


1 
72 


1 
2 


0.26 


0.11 


0.0064 


1.66 


0.28 



X,{0 = A'5-(C)-bA'5«(C)-f 0(m); 



(22) 
(23) 



namely, in the neighbourhood of jfi = the susceptibil- 
ity splits into a sum of two terms. The first of these 
expresses the contribution of the pion pole, 0(m~^), and 
can readily be expressed in a closed form 



2r.(C)^ 



mi 

m{C) 



(24) 



where the last equality is proved in Ref. [18]. The second 
term in Eq. (23), 0(to'^), is implicitly determined via 



i)S{q) 



(25) 



where the last line is obtained using {75, 7^} = 0. 

We can now proceed to our desired result. Equa- 
tion (25) entails 



<(C;m) = A'(0+O(m), 
where the vacuum chiral susceptibility is [9] 
d 



x{0 



dm{C) 



cr(m;C,A) 



(26) 



(27) 



Th=0 



Hence we arrive at a model-independent consequence of 
chiral symmetry and the pattern by which its broken in 
QCD; namely, 



^5(0 - 



m(C) 



^-(0 + 0(m) , 



(28) 



For illustration, in Table I we report numerical values 
computed from two models for the gap equation's kernel. 
Namely, we simplify the renormalisation-group-improved 
effective interaction in Ref. [22] 



Zig^Dp,{p-q)TUq,p) 



g{{p - qf) DlTiP - qy-rAq,p) , (29) 



wherein Dp^°{p—q) is the Landau-gauge free gauge-boson 
propagator, through the choice 



S ^6 



(30) 



which is a finite width representation of the form intro- 
duced in Ref. [23] . This interaction has been rendered as 
an integrable regularisation of 1/fc* [24]. Equation (30) 
delivers an ultraviolet finite model gap equation. Hence, 
the regularisation mass-scale can be removed to infinity 
and the renormalisation constants set equal to one. 

The kernel is completed by specifying the dressed- 
quark-gluon vertex. At leading-order in the systematic 
DSE truncation scheme [19, 20] the vertex is 



^a{q,p) 7a . 



(31) 



This defines the rainbow-ladder (RL) truncation. One 
can alternatively employ Ansdtze for the vertex whose 
diagrammatic content is unknown. A class of such mod- 
els, which has seen much use in diverse applications; e.g., 
Refs. [9, 21, 25, 26, 27, 28], can be characterised by [29] 



iT,{k, e) = ^SA(fc^ r) ja + {k + e)a 



^7-(fc + £)A^(fc2, 



where 



F{k^) - F{P) 

P --£2 



(32) 

(33) 
(34) 



with F = A,B; viz., the scalar functions in Eq. (3). This 
Ansatz satisfies the vector Ward-Takahashi identity and 
is often referred to as the BC vertex. 

Equation (28) is a remarkable result, which is nonethe- 
less readily understood. Recall that in the absence 
of a current-quark mass, the two- flavour action has a 
SUL(2)(g) SUfl(2) symmetry; and, moreover, that ascrib- 
ing scalar-isoscalar quantum numbers to the QCD vac- 
uum is a convention contingent upon the form of the 
current-quark mass term. 
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It follows that the massless action cannot distinguish 
between the continuum of sources specified by 

constant X J d'^x qix) e^'^-'^'^qix) , \e\ e [0,2n) . (35) 

Hence, the regular part of the vacuum susceptibility must 
be identical when measured as the response to any one 
of these sources, so that = X for all choices of 6. 
This is the content of the so-called "Mexican hat" po- 
tential, which is used in building models for QCD. The 
magnitude of X depends on whether the chiral symme- 
try is dynamically broken, or not; and the strength of 
the interaction as measured with respect to the critical 
value required for DCSB [9]. When the symmetry is dy- 
namically broken, then the Goldstone modes appear, by 
convention, in the pseudoscalar-isovector channel, and 
thus the pole contributions appear in X^ but not in the 
chiral susceptibility. It is valid to draw an analogy with 
the Weinberg sum rule [25, 30]. 

With Eq. (28) we have, in addition, provided a novel, 
model-independent perspective on a mismatch between 
the evaluation of the pion susceptibility using either a 
two-point or three-point sum rule. Namely, the two- 
point study [10] produces the pion pole contribution, X;^ , 



which is also the piece emphasised in Ref. [11], whereas 
a three-point method [8] isolates the regular piece, X^, 
because a vacuum saturation Ansatz is implemented in 
the derivation. Thus, the analyses are not essentially 
in conflict. Instead, they emphasise different, indepen- 
dent pieces of the susceptibility, which, with care, can be 
distinguished. However, in a sum-rules estimate of pion- 
nucleon coupling constants, only the regular piece should 
be retained [6]. 

We note in closing that the vacuum tensor susceptibili- 
ties, which can be related to the nucleon's tensor charges 
[31], can similarly be analysed. Such a study is underway. 
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